AD-A224  300 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-89-J-1530 
Task  No.  NR372-160 

TECHNICAL  REPORT  NO.  30 

FINITE  DEBYE- WALLER  FACTOR  FOR 
“CLASSICAL”  ATOM-SURFACE  SCATTERING 

by 

K.  Burke  and  W.  Kohn 

Prepared  for  publication 
in 

Phys.  Rev.  B  (1990) 

Department  of  Physics 
University  of  California,  Santa  Barbara 
Santa  Beirbcira,  CA  93106 


Approved  for  Public  Release 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States 
Government. 

This  document  has  been  approved  for  public  releeise  and  sale;  its  distribution  is  unlim¬ 
ited. 


August  1990 


REPORT  DOCUMENTATION  PAGE  ^  ^ 

read  INSTRUCT'.CjNS  ”1 

before  completing  form  I 

1.  HEKOH  1  NUMBER 

TECHNICAL  REPORT  NO.  30 

z.  govt  accession  no. 

i.  RECIPIENT’S  catalog  NUMBER  1 

N00014-01 

4.  title  (ond  Subllllo) 

FINITE  DEBYE-WALLER  FACTOR  FOR  "CLASSICAL" 
ATOM-SURFACE  SCATTERING 

s.  Type  of  report  .  period  covered 

TECHNICAL  REPORT 
1/01/90—12/31/90 

4  performing  oro.  report  num.EII 

T.  AUTNONf*; 

K.  BURKE  AND  W.  KOHN 

t.  CONTRACT  OR  grant  NUMBERlaJ 

N00014-89-J-1530 

>.  performing  ORGANIZATION  NAME  ANO  AOORCSS 

UNIVERSITY  OF  CALIFORNIA 

PHYSICS  DEPARTMENT,  SANTA  BARBARA,  CA  93106 
CONTRACTS  &  GRANTS,  CHEADLE  HALL,  ROOM  3227 

10.  program  element.  PROJECT.  TASK 

AREA  .  WORK  UNIT  NUMEERS 

TASK  NO.  372-160 

ELECTRONICS  &  SOLID  STATE  PHYSICS  PROGRAM 

800  N.  QUINCY,  ARLINGTON,  VA  22217 

IZ.  REPORT  DATE 

8/90 

IS.  NUMBER  OF  ^AGES 

14.  monitoring  AGENCY  NAME  4  kDORESS(ll  dllltront  Itom  Conirallln.  Olllcm) 

OFFICE  OF  NAVAL  RESEARCH  DETACHMENT 

1030  EAST  GREEN  STREET 

PASADENA,  CA  91106 

IS.  security  Class.  r>i  im*  r*PDr>j 

UNCLASSIFIED 

\Sm.  DECLASSIFICATION/  DOWNGBADING 
SCHEDULC 

fs.  distribution  statement  fot  thU  Report) 


"APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED" 


IT.  OlSTnieuTION  statement  (oI  ih»  Abtlracl  fifrtd  In  Block  30.  II  dllloioni  Irooi  Ktoorl) 

REPORTS  DISTRIBUTION  LIST  FOR  ONR 
PHYSICS  DIVISION  OFFICE— 

■  supplementahy  notes 


Phys.  Rev.  B  (.1990) 


If  KEY  WORDS  (Continue  on  #o»Ofoo  old#  ft  nocooo«rr  •’d  Identity  by  block  numbor) 

'''quar.Cur.  mechanical  scaCLerinc,  cia.-.si^ai  t;  r  c:  r  :  n  a  ,  cuasi-adiabat  ic  recimc, 
Debve-Wailer  Factor,  adiabatic  collision  variaoles,  atom-surtace  interaction, 
trajectory  approximation  „  ,, 


JO  abstract  /CwAllnuo  on  rovf/oo  ctJo  if  nocoooory  mnd  idontUr  by  block  numboti 

We  consider  atom-surface  scattering  at  low  surface  temperature,  fixed  incident 
energy  and  fixed  interaction  potential.  We  examine  the  limit  as  the  mass  m  of 
the  incident  particle  — »oon  and  thus  its  de  Broglie  wavelength  We  show 

that  in  this  -classical"  limit,  the  Debye- W^aller  factor  (a.  quantum  effect)  for  the 
strictly  elastic  scattering  amplitude  tends  to  a  finite  limit.  We  propose  scattering 
of  20  rneV  Ar  atoms  from  cold  Cu(lll)  as  a  promising  experiment  to  verify  this 
effect.  . 


1473  EDITION  Of  I  MOV  •»  1$  OSSOLETE 
5.'N  OlOJ-  LF.  OU-  6601 


nn 

ViJ  I  JAN  T1 


IECUPITY  CUASJIFICATION  OF  TMI$  PACE  Doim  tniorod) 


%  ■* 


Finite  Debye- Waller  Factor  for 
“Classical”  Atom-Surface  Scattering 

(Final  Version;  May  31,  1990) 

K.  Burke*  and  W.  Kohn 

Department  of  Physics 
University  of  California 
Santa  Barbara,  CA  93106 


Abstract:  We  consider  atom-surface  scattering  at  low  surface  temperature,  fixed 

incident  energy  and  fixed  interaction  potential.  We  examine  the  limit  as  the  mass  m  of 
the  incident  particle  -♦  oo,  and  thus  its  de  Broglie  wavelength  -►  0.  We  show  that  in 
this  “classical”  limit,  the  Debye- Waller  factor  (a  quantum  effect)  for  the  strictly  elastic 
scattering  amplitude  tends  to  a  finite  limit.  We  propose  scattering  of  20  meV  Ar  atoms 
from  cold  Cu(lll)  as  a  promising  experiment  to  verify  this  effect. 
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1.  Introduction 


Surface  scattering  experiments  are  generally  regarded  from  the  standpoints  of  two 
extreme  regimes:  the  quantum  mechanical  regime  or  the  classical  regime.  Consider  a 
surface  which  is  cold  (T  s:  0)  and  relatively  smooth  (exhibits  little  or  no  diffraction). 
Then  quantum  mechanical  scattering  is  exemplified  by  a  light  incident  atom  such  as  He 
with  a  low  initial  energy  of  say,  20  meV.  Under  these  conditions,  one-phodbn  peaks  may  be 
easily  distinguished,  and  the  technique  is  often  used  to  measure  surface  phonon  dispersion 
curves.^  On  the  other  hand,  classical  scattering  occurs  with  heavy  incident  particles  such 
as  Xe  at  hyperthermal  energies^  of  several  eV.  The  results  of  such  experiments  are  usually 
compared  to  either  classical  or  semiclassical  calculations.^  The  trsuosition  between  these 
two  regimes  is  continuous.  However,  it  is  generally  expected^  that,  as  the  mass  of  the 
incident  atom  is  increased,  characteristic  quantum  effects  disappear.  We  will  however 
show  that  (in  a  psurticular  sense)  quantum  effects  remain  finite,  if  the  imtial  energy  of  the 
incident  atom  is  kept  fixed. 

To  clarify  the  discussion,  we  re-express  the  distinction  of  regimes  in  terms  of  the  fol¬ 
lowing  dimensionless  parameters  :  A/d  and  n,  where  A  is  the  de  Broglie  wavelength  of  the 
incident  atom,  d  is  a  ‘‘range”  of  the  atom  surface  interaction,  and  fi  is  the  mean  number 
of  phonons  excited  by  the  collision.  Consider  first  scattering  from  a  rigid  lattice.  Then  the 
incident  atom  behaves  classically  if  A  C  d  and  quantum  mechanically  if  A  ^  d.  Similarily, 
the  lattice  behaviour  is  characterized  by  n.  If  after  the  collision  n  -w  0(1),  the  lattice 
exhibits  quantum  effects,  while  for  ft  ^  1,  it  behaves  classically.  For  the  experiments 
described  above,  both  incident  atom  and  lattice  behave  either  classically  or  quantum  me¬ 
chanically.  A  typical  energy  loss  spectrum  in  the  latter  regime,  integrated  over  all  outgoing 
directions,  is  shown  in  Fig.  1.  At  AE  ~  0,  there  is  a  substantial  no-loss  line  due  to  elastic 
scattering  from  the  surface  (this  has  been  broadened  in  the  figure  to  simulate  a  fimte  en- 
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ergy  resolution).  This  is  &  characteristic  quantum  effect  since,  for  classical  scattering  with 
a  finite  incident  energy,  there  is  necessarily  a  finite  energy  loss  to  the  lattice.®  We  refer  to 
the  scattering  as  weakly  inelastic  when  most  of  the  outgoing  distribution  lies  in  the  no-loss 
line. 

In  this  paper,  we  are  interested  in  a  mixed  regime  ;  A  C  d  but  fi  ~  0{1).  This 
may  be  achieved  physically  by  a  heavy  particle  (of  mass  m)  with  low  incident  energy 
E  striking  a  surface.  For  fixed  E,  as  m  is  increased,  the  initial  velodfy  of  the  incident 
atom  V  decreases  and  the  scattering  becomes  increasingly  adiabatic  and  more  and  more 
nearly  elastic.  In  the  limit  of  infinite  mass,  the  scattering  is  completely  elastic,  and  the 
adiabatic  limit  is  achieved.  For  large  but  finite  m,  the  scattering  is  almost  adiabatic,  and 
we  call  this  the  quasi-adiabatic  regime.  In  this  regime,  the  energy  loss  spectrum  of  Fig. 
1  is  compressed  toward  the  AE  =  0  region.  However,  assuming  we  have  sufficient  energy 
resolution  to  determine  its  shape,  there  are  three  possibilities  which  may  be  imagined  : 
(1)  the  scattering  becomes  completely  elastic  and  the  inelastic  shoulder  disappears,  (2)  the 
system  behaves  classically,  and  the  elastic  no-loss  line  disappears  or  (3)  a  finite  fraction  of 
both  elastic  and  inelastic  scattering  remain,  as  illustrated  in  Fig.  2.  We  show  in  this  work 
that  the  third  of  these  possibilities  actually  occurs,  providing  a  striking  demonstration 
of  a  quantum  mechanical  effect  that  does  not  vanish  as  the  mass  of  the  incident  atom 
grows  large.  (However,  the  splitting  between  the  elastic  peak  and  the  inelastic  shoulder 
approaches  zero.) 

In  section  II  of  this  paper,  we  make  a  realistic  estimate  of  the  Debye-Wjller  factor  in 
the  limit  of  large  m,  but  first  we  give  a  simple  plausibility  argument  for  the  above  claims. 
Let  us  denote  by  v  the  velocity  of  the  particle  during  the  interaction,  and  by  r(=  d/v)  the 
nominal  collision  time.  As  m  grows,  t  grows  as  The  slowly  varjring  force  exerted 

by  the  incident  atom  on  the  surface  atoms  can  excite  only  those  modes  with  frequencies 
u>  ^  1/r,  i.e.  only  long  wavelength  acoustic  phonons.  The  probability  of  exciting  a  single 
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phonon  of  frequency  u/  is  then  approximately^ 


~  Akluliu;), 


(1) 


where  is  the  perpen<iic\iiar  momentum  transfer  to  the  surface  for  elastic  scattering 

and  u|(a;)  is  the  contribution  from  a  normal  mode  of  frequency  u;  to  the  mean  square 
vibration  of  a  surface  atom  perpendicular  to  the  surface  in  the  undisturbed  crystal.  The 
elastic  fraction,  written  is  the  Debye-Waller  factor  .  For  small  3:  1  -  2W 

and  by  flux  conservation 

2W~.  Yi  (2) 

IV  —  1/r 


This  is  just  the  standard  expression  for  the  Debye- Waller  factor  (well  known  from  X-ray 
and  neutron  scattering  ^  ),  except  for  the  cut-off  at  ss  1/r  due  to  the  slowness  of  the 
incident  atom.  We  may  estimate  2W  using  the  bulk  properties  of  the  solid  in  the  Debye 
model:® 


2W 


{2kz) 


2ll  [ 
h  Jo 


2Mw  ) 


hu/JJ  )  ' 


(3) 


where  wjrj  is  the  Debeye  frequency  of  the  crystal,  M  the  mass  of  a  lattice  atom,  E_i  = 
k\/2m,  and  <  Uj  >  is  the  mean  square  of  zero  point  vibrations  in  the  crystal  perpendicular 
to  the  surface.  Note  that  this  is  independent  of  m,  and  therefore  W  remains  finite  as 
m  — »  oo.  Finally  we  note  that,  with  an  attractive  well  of  depth  D,  should  be  replaced 
by  E^-¥D  in  equation  (3)  (the  Beeby  correction).®  Thus,  we  still  find  a  finite  Debye- Waller 
factor  for  m  — »  oo. 


In  section  II  we  describe  our  model  for  the  atom-surface  interaction  and  the  approxi¬ 
mations  used  to  calcxUate  W  more  accurately.  The  quasi-adiabatic  regime  is  discussed  and 
illustrated  in  section  III  by  simple  models  for  the  interaction.  In  section  IV ,  we  apply  these 
results  to  noble  gas  scattering  from  metal  surfaces,  and  suggest  an  experiment  for  the  case 
of  Ar  scattering  from  Cu(lll).  We  choose  this  system  because  the  static  corrugation  is 
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i. 


expected  to  be  small  (so  that  diiTractioa  effects  may  be  neglected)  and  because  it  provides 
an  experimentally  accessible  system  with  a  relatively  heavy  incident  atom.  Properties  of 
other  systems  are  compared  with  this.  In  section  V  we  discuss  our  results  tmd  their  impli¬ 
cations  for  future  work,  both  theoretical  and  experimental.  Appendices  A  and  B  deal  with 
technical  details;  appendix  C  demonstrates  the  validity  of  the  trajectory  approximation 
for  weakly  inelastic  scattering  in  the  quasi-adiabatic  regime. 

II.  Theoretical  Formulation 

In  this  section,  we  introduce  a  model  for  the  atom-surface  interaction  and  demonstrate 
how  the  final  energy  and  angular  distributions  may  be  calculated  from  it.  We  are  interested 
in  the  quasi-adiabatic  regime,  achieved  by  m  — »  oo,  keeping  E  and  all  substrate  properties 
fixed.  However,  to  simplify  the  calctilations,  we  make  several  approximations  which  are 
valid  for  noble  gases  scattering  from  metal  surfaces  at  thermal  energies,  the  spediic  systems 
studied  in  this  paper.  For  example,  we  assume  throughout  that  the  displacements  of  the 
surface  atoms  are  small.  We  also  use  the  trajectory  approximation  (TA)  to  calc\ilate 
the  scattered  distribution.  Although  the  validity  of  the  TA  is  proved  here  (Appendix 
C)  only  for  weakly  inelastic  ,quasi-adiabatic  scattering  ,  we  expect  that  the  qualitative 
conclusions  remain  valid  also  for  strongly  inelastic  scattering.  Throughout,  we  use  the 
standard  notation  for  surface  scattering  :  A  vector  v  has  components  V  in  the  surface 
plane  and  vj  perpendicxilar  to  it;  its  length  is  denoted  by  v.  The  only  exception  to  this  is 
the  use  of  z  to  denote  the  component  of  the  incident  atom’s  position  perpendicular  to  the 
stirface.  We  also  use,  in  general,  primed  letters  to  denote  quantities  after  the  collision. 

A.  The  Atom-Surface  Interaction 

We  make  the  Bom-Oppenheimer  approximation^^  lor  the  electronic  motion,  which  is 
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valid  in  this  regime.  Then  a  general  atom-surface  interaction  may  be  written  V(r,  {u}) 
where  r  is  the  position  of  the  incident  atom,  and  the  set  {u}  is  the  set  of  all  displacements  of 
the  atoms  from  their  equilibrium  positions  in  the  crystal  lattice.  We  make  the  assumption 
(for  simplicity  only)  that  it  may  be  modeled  by  a  “perpendicular”  potential  function  of 
the  form 

V(r,{u})-VJz-Z(R,{u})),  (4) 

% 

where  r  =  (R,  2)  is  the  position  of  the  incident  atom  and  Z(R,  {u})  is  the  so  called  cor¬ 
rugation  function  for  the  atom-surface  system,  which  is  chosen  to  best  model  the  true 
many-body  potential.  The  corrugation  function  will  also  be  a  function  of  incident  energy, 
not  written  explicitly  here.  For  the  small  incident  energies  discussed  here,  the  incident 
atom  remains  far  from  the  surface  atoms,  and  this  approximation  will  be  valid.  For  small 
displacements  of  the  surface  atoms, 


V(r,{u})  =  V(r,(0})  ^  ^  u,.ViVj^ 


(5) 


where  1  labels  each  equilibrium  site  in  the  crystal,  and  {0}  indicates  no  displacements 
from  equilibrium;  the  gradient  with  respect  to  ui  is  denoted  Vj.  Then,  for  a  perpendicular 
potential,  as  defined  by  equation  (4), 


V(r,{u})  =  Vj.(j  -  Zo(R))  -  ^  ui.f(R  -  t) 

L 


dV 


dz 


*-Zo(R) 


+  . 


(6) 


where  Zo(R)(=  Z(r,  {0}))  is  the  static  corrugation  of  the  surface,  and  f(R)  is  the  so-called 
transfer  frmetion,  defined  by 


f(R)^ 


a^(R,{u}) 

duo 


(7) 


An  illustration  of  Zo(R)  and  f(R)  is  given  in  Fig.  3.  For  short-range  interactions,  dis¬ 
placements  of  atoms  beneath  the  surface  will  have  small  effects  on  the  potential  seen  by 
the  incident  atom  so  we  retain  only  the  Ij  =  0  (surface  layer)  term  in  the  sum  in  equation 
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(6).  From  equation  (6)  we  see  that  this  function,  coupled  with  knowledge  of  the  elastic 
scattering  potential,  "V ±{z  -  Zo(R)),  and  its  first  derivative,  is  sufficient  to  determine  the 
potential  for  small  displacements  of  surface  atoms.  Then,  for  a  system  of  negligible  static 
corrugation,  the  Hamiltonian  takes  the  form 


^  —  ^inc  +  5o  +  F^izUiR), 


(8) 


where  Rjne  =  p^/2m  -r  V’^(z)  is  just  the  Hamiltonian  for  the  incident  particle  scatter¬ 
ing  from  a  rigid  latttice,  Hq  is  the  Hamiltonian  for  the  undisturbed  crystal,  Fj_(z)  = 
—dVj_{z)l dz  is  the  force  exerted  on  the  incident  atom  by  the  crystal  and  ^(R)  =  53  ui,.f(R  —  L) 
is  linear  in  the  surface  atomic  displacements.  We  expand  the  displacements  in  terms  of 
the  harmonic  normal  modes  of  the  half-crystal: 


<? 


(9) 


where 

Ulu/n)  =  - ,  ,  ,=  .J=„:. 

In  this  expression,  advantage  has  been  taken  of  the  translational  symmetry  parallel  to  the 
surface  ;  N  is  the  number  of  atoms  in  the  crystal;  M  is  the  mass  of  a  crystal  atom;  Q  is  the 
phonon  momentum  parallel  to  the  surface;  q  =  {Q,93}  labels  the  3N  different  modes  and 
— g  =  {— Q,  —93};  <*g  creation  and  emnhilation  operators  for  phonons  labelled 


g,  with  frequency  u/g  and  displacement  vector  Vq  in  the  surface  layer.  The  function  u(u)g) 
is  just  the  magnitude  of  the  zero  point  displacement  in  the  bulk  due  to  a  bulk  phonon 
of  frequency  u>g,  while  surface  effects  are  contained  in  the  displacement  vector  Vg.  Its 
normalization  and  dependence  on  u/g  are  given  in  appendix  B.  We  also  define  the  surface 
Fourier  transform  of  the  transfer  function 
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(11) 


I 


where  A  is  the  area  of  the  surface  unit  cell,  and  S  represents  the  entire  surface.  The  sum 
over  L  now  yields 

{(R)  =  ^  r(G  ^  Q)  •  »uK)(a,  +  aL,).  (12) 

G.1 

As  only  low  frequency  modes  are  excited,  we  approximate  ?(K)  by  its  behaviour  in  the 
long  wavelength  limit  which  5rields 

«R)  =  +  a'.,)  ’  (13) 

■? 

(see  appendix  A). 

B.  The  Trajectory  Approximation 

Although  the  Hamiltonian  has  been  simplified  by  the  approximations  of  section  A, 
the  scattering  problem  is  still  not  easily  solved,  especially  for  strongly  inelastic  scattering. 
We  therefore  solve  it  within  the  trajectory  approximation  (TA),^°  in  which  the  incident 
particle  is  treated  classically,  while  the  surface  is  treated  quantum  mechanically.  For  weakly 
inelastic  scattering  in  the  quasi-adiabatic  regime,  this  is  shown  to  be  valid  in  appendix  C. 

First  we  calciilate  the  classical  trajectory  of  the  incident  atom,  with  the  atoms  fixed 
in  their  equilibrium  positions.  This  is  the  so-called  recoilless  classical  trajectory.  We  then 
have  a  simple  classical  one-dimensional  scattering  problem  with  a  potential  of  the  form 
shown  in  Fig.  4.  For  this  uniform  potential,  the  classical  equations  of  motion  decouple, 

and  we  write 

r  =  (V(,2rf(t))  (14) 

where  V  is  the  initial  velocity  parallel  to  the  surface  .  The  TA  Hamiltonian  then  becomes 
(see  equations  (8)  and  (13)) 

^TA  =  +  ^/2)  +  ^  "*■  “%)» 

q  q 

8 
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i 


where 

tr,(i)  =  (le, 


The  function  F(t)  is  just  evaluated  on  the  classical  trajectory  Zci(t).  This  Hamilto¬ 

nian  describes  a  set  of  uncoupled  linearly  driven  oscillators,  whose  final  states  and  corre¬ 
sponding  energy  and  momentum  changes  may  be  easily  calculated, The  distribution  of 
energy  loss  and  momentum  exchange  of  the  incident  atom  is  then  deduced  from  conserva¬ 
tion  laws  to  be 


-V'(A£:,AK) 


J  dt  J 


(17) 


where 

^/■(t.R)  =  e^cpj^  ;C;,{a;,)|7(exp[-iu»,f  -  iQ.R]  -  1)  (18) 

The  variables  AE  and  AK  denote  the  energy  loss  and  surface  momentum  exchange  re¬ 
spectively,  while  Uq{u})  is  the  Fourier  transform  of  Uq{t).  The  constant  term  (independent 
of  t  and  R)  in  equation  (18)  gives  a  contribution  to  the  final  spectrum  proportional  to 
S(AE)S(AK),  i.e.  it  is  the  Debye- Waller  factor  .  Using  equation  (16)  it  can  be  written  as 


2W  =  '£  +  Q.V)u“(u,,)|(v,)3|2.  (19) 

9 

This  result  has  the  same  form  as  equation  (3),  as  F'(O)  =  Akj  and  for  large  incident  mass 
F(<^q)  will  be  cut-off  at  about  1/r.  In  the  next  section  we  will  analyze  this  result  and  the 
energy  loss  spectrum  in  the  quasi-adiabatic  regime  in  more  detail. 


j 

in.  The  QuasUAdiabatic  Regime 


In  this  section  we  discuss  the  behaviour  of  the  energy  and  momentum  transfer  dis¬ 
tribution  found  in  section  II  (equation  (18))  in  the  quasi-adiabatic  regime.  We  introduce 
dimensionless  variables  to  show  explicitly  the  dependence  on  m  of  each  quantity  calcxilated. 
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In  partiCiilar,  we  find  that  the  mean  number  of  phonoas  excited  and  the  Debye-Waller  fac¬ 
tor  are  independent  of  m,  while  other  properties,  such  as  the  shape  of  the  energy  loss 
distribution,  scale  in  simple  ways  with  m. 

To  begin  with,  we  restrict  our  attention  to  normal  incidence  on  a  surface  at  zero 
temperature,  and  calculate  the  angle-integrated  energy  loss  spectrum.  For  the  classical 
trajectory  of  the  incident  atom,  we  choose  the  origin  in  space  and  time  to  be  at  the  position 
and  moment  of  turning  ,  i.e.  dzidt  =  0  at  z  =  0  and  <  =  0.  As  we  are  interested  in  the 
behaviour  of  the  system  as  m  -♦  oo  we  write  the  trajectory  in  terms  of  the  dimensionless 
variables  used  in  section  I:(^  =  z/d,  j=t/r  and^2  ^  v^{  z)/E.  The  equation  of 
motion  for  ^  is  then 

(20) 

where  ("  =  d^^/ds^.  The  boundary  conditions  are  C^(±oo)  =  ±1.  With  this  choice  of 
origin,  ^(s)  is  symmetric  in  s,  and  independent  of  m  (but  not  of  E).  In  terms  of  these 
variables,  the  energy  loss  spectrum  is  denoted  where  e  =  ^Et /h  and  the  normalization 
is  chosen  so  that  /  A'<d«/2ir  =  1.  Then,  from  equation  (16),  with  R  =  0,  we  may  write  its 
Fourier  trimsform  as 


^/{s)  =  exp{a[x(s)  -  x(0)l}. 

(21) 

x(»)= 

Jo 

(22) 

„  /  E  V  <  >, 

(23) 

and  Ci,'  is  the  Fourier  transform  of  C"(s)-  In  these  expressions,  we  use  the  long  wavelength 
limit  to  calc\ilate  the  phonon  properties,  as  described  in  appendix  B.  The  quantity  <  Uj  >, 
is  a  measure  of  the  zero  point  displacements  perpendicular  to  the  surface  of  atoms  in  the 
surface,  and  is  also  described  in  appendix  B.  The  Debye- Waller  factor  is  then  given  by 
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exp-ax(O),  so  that^^ 


where  In  =  This  shows  explicitly  that  2W  is  independent  of  m,  and  has 

the  form  of  equation  (3). 


To  illustrate  these  results  we  use  two  simple  potentials  for  which  analytic  results  may 
be  obtained  :  the  exponential  repulsion  and  the  Morse  potential.  For  both  potentials  we 
choose  d  to  be  the  decay  length  of  the  repulsive  part  of  the  potential  and  a  =  0.1,  a 
value  appropriate  for  the  experimental  conditions  discussed  in  section  IV.  For  the  Morse 
potential,  we  also  choose  a  well-depth  D  =  2.ZE  for  the  same  reason.  The  dimensionless 
forms  for  the  Fourier  transform  of  the  acceleration  are  then 


2irv 

sinhTTi/ 


(25) 


for  the  pure  exponential  repulsion  and 

_  2iru  cosh/?7ri/ 
sinhTTi/  coshTTt/ 

for  the  Morse  potential.  Here  is  a  function  of  E/ D,  with  a  value  of  1.6  for  our  case.^* 
These  frmctions  are  plotted  in  Fig.  4,  and  their  corresponding  final  energy  spectra  in  Fig. 
5.  The  dramatic  difference  is  due  to  the  attractive  well,  which  accelerates  the  incident 
particle  as  it  passes  over  it.  From  Table  I,  we  see  that  this  greatly  enhances  Ii,  leading 
to  a  major  decrease  in  the  Oebye- Waller  factor.  We  stress  the  fact  that  this  reduction  is 
entirely  due  to  the  attractive  well,  and  is  not  due  to  the  large  incident  mass. 


To  clauify  further  our  discussion  of  the  loss  spectrum,  we  introduce  a  number  of  pa¬ 
rameters  characteristic  of  the  collision.  From  the  spectrum  given  by  equation  (18),  it 
can  be  shown  that  the  mean  phonon  frequency  to  which  the  incident  atom  couples  is 
jj  =  {I2I I\)'^~^ •  The  mean  energy  loss  may  also  be  calculated,  and  we  find  ^E  =  haJ2/T. 
Then  a  measure  of  the  number  of  phonons  excited  during  the  collision  is  given  by  ft  = 
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ilEfuy  or 


n  =  2W.  (27) 

The  meaning  of  this  result  is  intuitively  clear  :  when  few  phonons  me  excited  (n  <  1), 
the  scattering  is  predominantly  elastic  (weakly  inelastic  scattering  ),  while  for  ft  >  1,  the 
converse  is  true. 


Next,  we  briefly  discuss  temperature  eflects.  Equation  (18)  may  be  easily  generalized^*^ 
to  finite  s\irface  temperature,  T.  The  only  change  this  leads  to  in  equation  (19)  is  an 
extra  mviltiplicative  factor  of  coth.(ujq/2kBT),  due  to  the  thermal  motion  of  the  crystal 


atoms.  For  small  T,  the  Debye-Waller  factor  behaves  as  a  gaussian  function  of  T.  We  find 
w(T)  =  w{Q)  +  with 


where  k^Tr  =  h/r.  From  the  definition  of  a  (equation  (24))  amd  using  r  =  d/v,  it  is 
straightforward  to  show  that  T  is  independent  of  the  functional  form  of  the  potential. 
However,  the  range  of  temperature  over  which  the  Debye- Waller  factor  displays  gaussian 
behaviour  is  determined  by  the  form  of  the  interaction  potential.  For  high  temperatures 
{T  Tt)  the  usual  exponential  dependence  of  the  Debye- Waller  factor  is  regained,  with 


a  decay  temperature  of 


(29) 


The  dependence  of  T  on  the  other  parameters  differs  markedly  from  standard  expressions,^ 
due  to  the  cutoff  in  Fx(w)  at  1/r. 


For  any  finite  temperature,  as  long  as  m  is  finite,  there  is  a  finite  Debye- Waller  factor 
.  However,  as  m  ->  oo  with  T  finite,  T  »  Tr,  so  that  coth(u>g/2ifej3r)  >  1  for  all  excited 
phonons,  and  the  probability  of  excitation  is  greatly  enhanced.  Hence,  as  m  — »  oo,  the 
Debye- Waller  factor  approaches  zero. 
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The  angular  distribution  is  gaussian,  of  width 


^  /8a/3/i,  Va , 

(30) 

where  u®  =  fi/md  and  ct  is  the  transverse  speed  of  sound  in  the  crystal.  The  dimensionless 
constant  /i  is  an  eljistic  property  of  the  crystal,  with  a  value  of  0.9  for  Cu  (Appendix  B). 


rV.  Application  to  Scattering  of  Rare  Gases 


We  begin  our  discussion  of  results  with  the  system  ;  Ar  incident  on  Cu(lll)  with  an 
initial  energy  of  20  meV.  This  is  a  good  system  for  experimental  verification  of  this  effect 
for  a  number  of  reasons.  Of  primary  importance  is  the  fact  that  Ar  is  stifficiently  heavy,  so 
that  the  quasi-adiabatic  regime  may  be  achieved  at  low  incident  energy.  Furthermore,  the 
well-depth  of  the  interaction  potential  is  sufficiently  small  so  that  the  elastic  peak  remains 
resonably  strong  (see  Fig.  4).  Another  useful  feature  which  simplifies  the  calculation  is  the 
smoothness  of  a  (111)  metal  surface,  which  allows  us  to  neglect  diffraction.  However,  even 
with  these  simplifications,  we  offer  only  a  semiquantitative  calculation  of  the  Debye- Waller 
factor  .  Our  purpose  is  simply  to  show  that  this  effect  occurs  for  Ar  scattering  &om  Cu, 
and  that  it  is  detectable. 


For  the  elastic  scattering  potential  we  asstime  a  simple  exponential  repulsion  plus  a 
long  range  Van  der  Waals  attraction: 

(31) 


The  constants  C3  and  ^ref  given  by  Zaremba  and  Kohn.^®  We  expect  the  Ar  -  Cu 
repulsion  to  have  approximately  the  same  z-dependence  as  the  He  -  Cu  repulsion  (as 
suggested  by  the  Norskov  recipe)^®.  The  proportionality  constant^^  is  4.6.  We  also  take 
the  He  -  Cu  repulsion  from  Zaremba  and  Kohn.^®  The  values  we  use  are  listed  in  Table 
II,  along  with  other  characteristics  of  the  potential.  These  numbers  are  only  estimates, 
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and  shoxild  not  be  regarded  as  highly  accurate.  The  most  important  characteristics  of  the 
potential  are  the  decay  length  of  the  repulsion  and  the  well-depth,  which  is  45  meV.  For 
our  purposes,  its  properties  are  very  similar  to  that  of  a  Morse  potential  of  equal  depth 
and  decay  length,  the  difference  in  the  long-range  attraction  being  of  no  importance. 

The  energy  loss  spectrum,  integrated  over  all  outgoing  directions,  calculated  in  the  TA, 
is  plotted  in  Fig.  5.  The  surface  is  at  zero  temperature.  The  spectrum  has  been  convoluted 
with  a  gaussian  of  width  2  meV,  to  simulate  an  experimental  resolution  of  10%.  Its  shape 
is  very  similar  to  that  of  the  Morse  potential  shown  in  Fig.  4,  except  for  the  cutoff  at 
the  incident  energy.  The  elastic  peak  contains  about  10%  of  the  scattered  particles,  while 
another  50%  are  scattered  into  the  inelastic  shoulder.  A  further  40%  remain  stuck  to  the 
surface  (negative  final  energy  in  the  trajectory  approximation).  The  energy  of  the  incident 
particles  in  the  interaction  well  is  ss  65  meV,  whereas  the  calculated  mean  energy  of  the 
scattered  particles  is  ss  46  meV,  fortunately  in  qualitative  agreement  with  the  recoilless  TA, 
for  which  both  energies  would  be  65  mcV.  (Qualitative  considerations  lead  us  to  believe 
that  a  better  calcxilation  of  the  Debye- Waller  factor  would  yield  a  fraction  somewhat  larger 
than  the  10%  fraction  quoted  above.) 

We  conclude  our  discussion  of  Ar  on  Cu  by  examining  other  properties  of  the  outgoing 
distribution.  From  Table  III,  we  see  that  the  angular  distribution  is  q\iite  narrow  (a  half¬ 
width  of  2®),  so  that  a  detector  with  a  wide  aperture  (about  8®)  could  count  all  scattered 
particles  simultaneously.  We  can  also  show  that  the  scattered  distribution  is  insensitive 
to  the  angle  of  incidence,  so  the  experiment  could  be  performed  with  the  convenient  total 
scattering  angle  of  90°.  Lastly,  the  temperature  dependence  of  the  no-loss  peak  of  Fig.  5 
is  given  in  Fig.  6.  As  noted  in  section  III,  it  is  a  gaussian  near  T  =  0,  but  becomes  an 
exponential  for  higher  temperatures.  The  figure  shows  that  the  surface  needs  to  be  cooled 
to  a  temperature  T  S  SOif  to  ensure  that  the  elastic  peak  is  measurable. 

We  have  also  examined  the  possibility  of  using  other  heavy  noble  gases  as  incident 
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atoms.  Clearly,  for  Kr  or  Xe,  the  scattering  is  closer  to  adiabatic  than  for  At.  However,  the 
larger  well-depths  which  these  incident  atoms  experience  reduces  the  Debye- Waller  factor 
so  much  that  the  no-loss  line  is  very  small.  We  have  also  examined  different  crystals,  to 
see  if  a  larger  Debye- Waller  factor  is  likely  to  be  found.  In  Table  IV,  we  give  results  for 
three  other  surfaces  ;  Ag,  A1  and  Au.  The  potentials  and  surface  elastic  properties  for 
these  surfaces  were  estimated  in  the  same  way  as  those  of  Cu.  Although  the  well-depth 
for  Ag  is  very  close  to  that  of  Cu,  the  increase  in  a,  due  to  the  lower  Debye  temperature, 
reduces  its  Debye- Waller  factor  .  In  the  case  of  Al,  it  is  the  increased  well-depth  that  leads 
to  the  reduction  in  the  small  Debye- Waller  factor  ,  while  for  Au,  both  these  effects  occur. 
However  any  of  these  surfaces  are  possible  candidates  for  seeing  this  effect.  On  the  other 
hand,  the  well-depth  for  Ar  on  W  is  so  leirge^®  (—  100  meV)  that  it  is  an  unpromising 
choice. 


V.  Conclusions 

A  heavy  incident  atom,  with  fixed  (thermal)  energy  scatters  almost  adiabatically  from 
a  target  surface.  For  T  =  0,  in  the  limit  of  infinite  mass,  the  Debye- Waller  factor  tends  to 
a  finite  limit  (<  1).  This  is  because  the  slowly  moving  incident  atom  can  only  excite  low 
frequency  phonons.  We  used  the  trajectory  approximation  to  estimate  the  final  energy 
spectrum.  We  find  that  20  meV  Ar  striking  a  cold  (T  ^  30K)  Cu(lll)  surface  should 
exhibit  substantial  (about  10%)  purely  elastic  scattering.  Other  rare  gas/metal  systems 
appear  to  be  less  favorable.  This  Debye- Waller  factor  becomes  negligible  at  T  s:  100  K. 

We  have  also  shown  in  Appendix  C  that,  for  quasi- adiabatic  weakly  inelastic  scattering, 
the  trajectory  approximation  is  valid.  The  question  of  its  validity  for  strongly  inelastic 
scattering  and  for  the  scattering  of  light  atoms  will  form  the  subject  of  future  work. 

Although  we  used  several  simplifying  assumptions  we  stongly  believe  that  our  basic 
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conclusion,  that  in  the  adiabatic  limit  (m  — *  oo)  the  zero  temperature  Debye- Waller  factor 
tends  to  a  finite  value,  is  independent  of  these. 
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Appendices 


A.  The  Surface- Averaged  Transfer  Function. 


In  this  appendix,  we  derive  general  expressions  for  the  surfaced-averaged  transfer  func¬ 
tion  for  a  given  atom-surface  interaction.  We  restrict  ourselves  to  potentials  which  may 
be  written  as  a  sum  of  pair  potentials,  and  whose  static  corrugation  is  relatively  smji.ll 
We  also  consider  only  monatomic  surfaces.  We  find  the  perpendiculaf  component  is  a 
constant,  independent  of  the  interaction,  and  the  parallel  components  vanish. 

For  a  many-body  interaction  which  is  the  sum  of  pair  potentials,  we  may  write 


^.(;r-n-u,|)^Vi(»-Z(R,{u}))  (A.1) 

1 

where  Vj_(z  —  Z(R,  {u}))  is  the  perpendicular  potential  approximating  the  true  potential 
and  v(r)  is  the  interaction  between  the  incident  atom  and  a  single  surface  atom,  a  distance 
r  apart.  The  label  1  ranges  over  all  equilibrium  positions  F]  in  the  crystal,  while  U]  are 
the  displacements  of  the  crystal  atoms  from  these  positions.  Now  consider  an  infinitesimal 
displacement  u  of  the  1  =  0  atom  (on  the  surface  ).  Then 


VJz  -  ZoiR)  -  u.f(R))  =  ^  v(!r  -  n  -  «l,o.u|) 

I 


(A.2) 


where  f(R)  is  the  transfer  function  (defined  in  the  main  text,  equation  (7)).  Expanding 
the  potential  functions  about  their  equilibrium  values,  we  find 


f(R)  = 


v'(r) 


(A.3) 


V'^(z  -  Zo(R)) 

where  r  =  r/r.  The  primes  on  the  potentials  indicate  the  derivatives  with  respect  to  their 
explicit  dependent  variable  ,  i.e.  v'(r)  =  dv/dr,V'^(z  -  Zq(R))  =  dV^/dz.  Note  that, 
despite  appearances,  f(R)  is  implicitly  independent  of  z,  by  definition.  Then 


v'(r) 


^o(R)) 


(A.4) 
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Now  we  make  use  of  our  second  assumption,  viz.  that  Zo(R)  is  fairly  smooth,  to  replace 
V'j_(z  -  Zo(R))  in  the  above  expression  by  its  average  and  to  bring  it  outside  the  integral 


f(0)  =  f  (A  5) 

AV'  Js 

where 

=  jJ^d^RV'jz  -  Zo(R))  (A.6) 

the  average  being  taken  over  only  one  unit  cell  as  the  static  potential  has'the  translational 

symmetry  of  the  surface  .  From  equation  (A. 4)  we  see  from  symmetry  that  the  transverse 
component, 

f{0)=0.  (A.7) 

To  find  /3(0),  we  note  from  equation  (A.6) 

=  (A.8) 

using  equation  (A. 2)  with  u  =  0.  By  replacing  r  -  ril  with  a  new  variable  r  in  each  sum 
we  find 

=  4  /  (A.9) 

-  A  Js  r 

which,  when  substituted  into  equation  (A. 5),  yields 

/3(0)  =  1  (A.IO) 


the  desired  result. 

B.  The  surface  displacement  partial  spectral  density  perpendicular  to  the  surface  for  an 
isotropic  elastic  medium. 

the  adiabatic  limit  the  incident  atom  excites  only  the  low  frequency  phonons  of 
the  crystal,  so  we  use  only  the  long  wavelength  properties  of  the  crystal.  We  estimate  th' 
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desired  quantities  by  treating  the  crystal  as  isotropic,  which  considerably  simplifies  the 
calculations. 


The  properties  of  an  isotropic  elastic  medium  with  a  surface  (at  z  =  0)  were  originally 
caJculated  by  Maradudin  and  Mills.  We  define  to  be  an  eigenvector  of  the  dynamical 
matrix  of  the  (half-)crystal,  with  normalization 

1 

where  1  labels  the  equilibrium  positions  of  the  atoms  in  the  crystal,  and  N  is  the  number 
of  atoms  in  the  crystal.  As  the  crystal  has  translational  symmetry  parallel  to  the  surface, 
we  may  write 

v,(l)  =  v,(/3)e*QI'.  (B.2) 

For  bulk  phonons  far  within  the  crystal  — »  oo),  =  e^e^*^*,  where  is  a  polar* 

ization  vector  of  unit  length.  For  btdk  phonons  at  the  surface  of  the  crystal  (/j  =  0),  the 
lengths  of  these  vectors  are  greater  than  1,  due  to  enhancement  of  atomic  displacements 
at  the  surface.  Finally,  for  surface  phonons,  it  should  be  noted  that,  for  low  frequency 
displacements,  as  u;  — ♦  0,  the  displacement  vectors  grow  as  (for  bulk  phonons,  they 

MC  independent  of  ci/).  In  terms  of  these  displacement  vectors,  we  may  now  define  the 


Green’s  function  by 


-  u;| 


(8.3) 


which  has  been  Fourier  transformed  parallel  to  the  surface.  Here  a  and  /3  denote  the 
Cartesian  components.  Then  the  displacement  partial  spectral  density  perpendicular  to 


the  surface  is  defined  as 


yielding 


£  d^Q  2u» 

iV  Jq  (25r)2  TT 


Imd33(Q,l3,/3;u/  -  ie) 


(B.4) 


(B.5) 
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where  Kq  is  the  two-dimensional  Debye  wavenumber,  defined  by 


where  p  is  a  dimensionless  function,  given  (for  =  0)  in  equation  (A6)  of  the  appendix  of 
Stutki  and  Brenig,^*^  and  is  the  transverse  velocity  of  sound  in  the  bulk.  Substituting 
this  expression  back  into  equation  (B.S),  and  using  the  definition  of  Kq  given  in  equation 
(B.6)  and  the  standard  Debye  frequency,  we  find 


where 


and 


i;  w)  = - = — -  w) 


lih  =  0;u.)  =  — ^:g<“>(oo)  - 

i  -r  (X 


=  f 

Jo 


pih  = 

~i2+2n 


(B.8) 

(B.9) 

(B.IO) 


Here  <r  =  ct/ci  is  the  ratio  of  transverse  and  longitudintd  velocities  in  the  bulk.  Further¬ 
more,  as  the  Rayleigh  velocity  is  the  minimum  velocity  for  any  elastic  wave  in  the  surface, 
G(*)  =  0  for  X  <  cjilet,  where  cji  is  the  Rayleigh  wave  velocity.  Thus  for  low  frequencies 
=  7(^3)  •  constant.  From  equation  (B.S)  it  is  clear  that  7(13)  — *  1  as  13  — ►  00  (the 

interior  of  the  crystal),  while  7,  =  7(^3  =  0)  is  generally  larger  than  one,  due  to  enhance¬ 
ment  of  atomic  displacements  at  the  surface.  Furthermore  we  write  <  Uj  >,=  7*/(4Afu»£j) 
as  a  measure  of  the  mean-square  displacement  perpendicular  to  the  surface  at  the  surface 
(Af  is  the  mass  of  a  surface  atom).  These  definitions,  and  the  spectral  density  of  equation 
(B.S),  lead  to  the  simple  form  of  equations  (22)  and  (23)  of  the  main  text. 

20 


I 


The  result  that  5*(u;)  is  proportioaal  to  for  small  u  clearly  depends  only  on  the  three- 
dimensional  nature  of  the  crystal.  It  is  only  in  estimating  the  proportionality  constant 
that  we  use  the  isotropic  model  as  an  approximation.  For  anisotropic  crystals,  we  write 
a  =  Ci/di,  where  the  velocities  have  been  averaged  over  all  directions  in  the  bulk.^^ 

We  conclude  by  noting  that  the  factor  fj,  occuring  in  equation  (30)  of  the  text,  and 
contributing  to  the  width  of  the  angixlar  distribution,  is  simply  C»(^)(oo  oo). 

C.  Validity  of  the  Trajectory  Approximation  in  the  Quasi-Adiabatic  Regime 

In  this  appendix  we  show  that,  for  weakly  inelastic  scattering  and  small  displacements 
of  the  surface  atoms,  the  recoilless  TA  is  valid  in  the  quasi-adiabatic  regime,  i.e.  when 
the  mass  of  the  scattered  particle  becomes  large.  This  is  because  the  incident  atom  does 
indeed  behave  classically,  although  the  surface  does  not.  We  calculate  the  angle-integrated 
energy  loss  spectrum  using  the  distorted  wave  Bom  approximation  ,  and  compare  resxxlts 
with  those  foxmd  from  equation  (21). 

In  the  language  of  quantum  mechanical  scattering  theory  we  write,  for  the  energy  loss 
spectrum,^^ 

^(<)  =  E  E I  <  “i A'-kio  >!’<(£»  -  Co  - 1)  (c.i) 

where  is  the  T  matrix  for  scattering  the  incident  atom  with  wavevector  k  to  final 

wavevector  k\  with  energy  loss  «.  It  is  an  operator  on  the  states  ja  >  of  the  undisturbed 
crystal,  where  [0  >  is  the  ground  state.  These  states  have  energy  Ea-  To  find  this 
matrix  element,  we  linearize  the  potential  in  the  displacements  of  the  surface  atoms  (see 
equation  (6))  and,  for  weakly  inelastic  scattering,  apply  the  first  order  distorted  wave  Bom 
approximation:^^ 

<l 
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where  is  the  leading  inelastic  contribution  to  the  T  matrix.  Here  ^4,(2)  is  the  scat¬ 
tering  solution  for  an  incident  atom  of  perpendicular  wavevector  component  ^3  scattered 
from  V()(«).  We  have  again  used  r(G  -r  Q)  =  (0,^g,o)  simplicity.^^ 

The  delta  functions  in  equations  (C.l)  and  (C.2)  impose  conservation  of  energy  and 
surface  crystal  momentum  on  the  outgoing  states.  At  fixed  incident  energy,  as  m  — »  00, 
both  and  — *  00,  so  that  WKB  wavefunctions  may  be  used  to  calctilate  the  matrix 
element  containing  We  write  * 

<  >=  r  r d^nd^)FAz)^h,{^)-  (C.3) 

*  J-ao  *  Jz,  * 

Here,  the  origin  of  the  z-axis  is  chosen  at  the  classical  turning  point  of  the  particle  with 
energy  E,  while  z^  is  chosen  such  that  z©  ^  A,  but  also  so  that  z©  — »  0  as  m  — »  00  (in 
a  fashion  described  below).  Thus,  to  evaluate  the  second  integral,  we  may  use  the  WKB 
states.  We  define  these  by 

where 

<^(2)  =  r  fc3(z')dz'  (C.5) 

and  fc3(i)  =  ^k\  -  2mViX^).  Then  the  physical  scattering  wavefunctions  may  be  written 

♦*.(*) = (c-6) 


The  second  integral  appearing  in  equation  (C.3)  will  have  four  contributions.  Consider  the 
first  of  these  : 


We  may  expand  Al;3(z)  (=  ^3(2)  —  ^3(2))  fof  l^rge  m,  and  find  Afc3(z)  =  mulki^z)  to 
leading  order  in  m^/^.  In  this  expression,  the  denominator  never  vanishes,  as  we  have 
excluded  the  turning  point  from  our  integration  interval.  Next  we  consider  the  classical 
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trajectory  Zcl(0  associated  with  the  static  perpendicular  potential.  We  define  it  as  a  time- 
symmetric  function  starting  at  -oo,  reaching  the  turning  point  at  t  =  0  and  continuing 
on  to  -1-00.  On  the  interval  0  to  oo,  it  is  an  invertible  function  of  t,  so  that  t  — 
is  well-defined.  Changing  variables  from  z  to  t  using  this  function,  we  find  a  very  simple 
result  for  the  phase  of  our  integral,  viz. 


(C.8) 

% 

where  to  =  same  change  of  variables  for  the  entire  integral,  and  use 

jfe3  ~  jfe3  in  all  other  functions  to  get  the  leading  order  contribution: 

(C.9) 

where  v  is  the  initial  velocity  of  the  incident  particle.  Finally,  as  m  — »  oo  and  A  — *  0,  so 
that  to  — >  0,  the  integrand  is  well-behaved  and  we  have 

(C.IO) 

where  F^(u;)  is  the  Fourier  transform  of  the  classical  force  exerted  on  the  incident  atom.  In 
an  obvious  notation,  we  find  =  0  due  to  the  rapid  oscillation  of  the  phases 

in  the  integrals,  while,  from  symmetry,  Af^““)  =  |i4^  ^l^vFj^(u;)/2.  Then,  conservation 
requirements  yield 

<  >=  (C.ll) 

Inserting  thus  matrix  element  back  into  the  T  matrix,  and  using  that  to  calcxilate  the 
energy  loss  spectrum,  we  find  identical  results  to  the  TA  for  weakly  inelastic  scattering 
(i.e.  expanding  the  exponential  in  equation  (18)  in  a  power  series,  and  keeping  only  the 
first  term). 

We  conclude  by  showing  that  the  first  integral  in  equation  (C.3)  is  negligible  in  the 
adiabatic  limit.  First  note  that,  for  z  <  —zq,  we  may  again  use  WKB  to  show  that  this 
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contribution  becomes  exponentially  small  as  m  — »  oo.  Then,  around  the  turning  point 
z  =  0,  there  is  a  small  interval,  -d  to  <f,  in  which  the  potential  may  be  well-approximated 
by  a  linear  function.  The  value  of  d  depends  only  on  E  and  the  shape  of  the  potential,  and 
is  therefore  independent  of  m.  For  -d  <  z  <  d,  we  may  use  Airy  functions  to  approximate 
the  wavefunctions.  If  we  choose  Zq  =  {2mF_{t  =  0))“^'^^,  where  F_(0)  is  the  classical  force 
at  the  turning  point,  we  find  the  integral  from  -Zq  to  Zq  to  be  finite,  and  proportional  to 
Zq.  This  choice  obeys  the  requirements  stated  above,  and  )delds  a  vanishing  contribution 
(of  0(m~^/^))  as  m  — ♦  oo.  This  completes  the  proof. 
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Table  Captions 


(I)  Constants  characterizing  quasi-adiabatic  scattering  for  three  different  potentials: 
the  exponential  repulsion,  the  Morse  potential  and  the  potential  used  to  model  the  Ar/Cu 
system. 

(II)  Properties  of  Cu  surface  and  Ar-Cu  interaction  used  to  calcufate  Fig.  5.  The 
sources  are  given  in  the  main  text.  The  potential  parameters  should  not  be  regarded  as 
highly  accurate. 

(III)  Cedciilated  properties  of  quasi-adiabatic  scattering  for  the  energy  loss  spectrum 
of  Fig.  5. 

(IV)  Calculated  Debye-Waller  factors  for  Ar  on  three  different  crystals.  D  denotes  the 
estimated  weli-depth  and  a  is  defined  by  equation  (23). 
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constant 

exponential 

Morse 

Ar/Cu 

1 

h 

2.07 

25.1 

21.7 

h 

0.72 

31.0 

23.0 

h 

0.41 

52.5 

34.4 

h 

0.31 

111 

64.9 

^-2W 

0.93 

0.04 

0.10 

fl 

0.07 

3.1 

2.3 

Table  I 


parameter 

symbol 

strength  of  repulsion 

85.7  eV 

repulsive  decay  length 

d 

0.37A 

Van  der  Waals  constant 

Ci 

1.50  eV  A’ 

position  of  reference  plane 

Zrtf 

0.259  A 

depth  of  attractive  well 

D 

45.4  meV 

surface  phonon  enhancement 

7 

1.8 

Table  H 
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1 

metal  | 

1 

D(meV) 

c»/aCu 

<% 

1 

■  '"I 

Ag 

44.6 

1.60 

3.2 

A1 

64.5 

1.18 

1.1 

Au 

69.6 

2.17 

0.09 

Table  IV 
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Figure  Captions 


(1)  Typical  (angle-integrated)  energy  loss  spectrum  for  atom-scattering 
from  a  cold  surface.  The  no-loss  line  heis  been  broadened  to  simulate  finite 
energy  resolution. 

(2)  Change  in  energy  loss  spectrum  for  a  heavy  particle  when  incident 
mass  is  doubled.  The  inelastic  shoulder  grows  in  height  (by  a  factor  of  n/2) 
and  is  compressed  towards  A£  =  0  (again  by  a  factor  of  so  that  its  total 
area  remains  unchanged.  The  no-loss  peak  does  not  change  at  all,  so  that 
the  Debye-Waller  factor  remains  finite. 

(3)  Schematic  of  corrugation  function  for  a  many-body  potential,  (a) 
The  heavy  lines  indicate  atomic  positions  and  corrugation  function  at  equi¬ 
librium  positions.  The  dashed  line  shows  the  change  in  the  corrugation 
function  due  to  moving  the  central  atom  a  small  amount  perpendicular  to 
the  surface,  (b)  The  transfer  function  in  the  perpendicular  direction. 

(4)  Fourier  transform  of  classical  accelerations  for  two  potentials  in  di¬ 
mensionless  variables  :  u  =  wr  and  C”  =  iwh. 

(5)  Energy  loss  spectra  for  two  potentials  in  the  quasi- adiabatic  regime. 
Only  weakly  inelastic  scattering  occurs  for  the  exponential  (n  =  0.07),  but 
for  the  Morse  potential,  the  acceleration  due  to  the  well  leads  to  strongly 
inelastic  scattering  (n  =  3.i). 

(6)  Calculated  total  final  energy  spectrum  for  20  meV  Ar  striking  Cu 
(111)  at  zero  temperature.  The  scattering  is  about  10%  elastic.  An  energy 
resolution  of  only  10%  was  assumed,  and  the  no-loss  peak  is  clearly  visible. 

(7)  The  dependence  on  surface  temperature  of  the  peak  shown  in  Fig. 
5. 
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